Toward scalable information processing with ultracold polar molecules in an electric 

field: a numerical investigation 
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We numerically investigate the possibilities of driving quantum algorithms with laser pulses in a 
register of ultracold NaCs polar molecules in a static electric field. We focuse on the possibilities of 
performing scalable logical operations by considering circuits that involve intermolecular gates (im- 
plemented on adjacent interacting molecules) to enable the transfer of information from one molecule 
to another during conditional laser-driven population inversions. We study the implementation of 
an arithmetic operation (the addition of or 1 on a binary digit and a carry in) which requires 
population inversions only and the Deutsch-Josza algorithm which requires a control of the phases. 
Under typical experimental conditions, our simulations show that high fidelity logical operations 
involving several qubits can be performed in a time scale of a few hundred of microseconds, opening 
promising perspectives for the manipulation of a large number of qubits in these systems. 



I. INTRODUCTION 

The past few years have witnessed remarkable exper- 
imental achievements on the realization of elementary 
quantum logical operations on physical systems [lHSl- 
Building a quantum processor that could handle a large 
number of quantum bits (qubits) represents the next 
milestone to reach toward the realization of a practical 
quantum computer, but no technology is currently avail- 
able for processing about one hundred of qubits which 
is the expected minimum number of qubits required to 
overcome powerfull current classical computers. Due to 
their rich inner energy structure that can be used to en- 
code information, molecules offer promising propects for 
scalable quantum information processing and have at- 
tracted a lot of attention recently. Followingthe pionner- 
ing work of de Vivie-Riedle and coworkers [1, 0], several 
groups have explored the possibily of encoding qubits in 
ro- vibrational states of a single diatomic molecule [8l4l8j 
or polyatomic molecul e I19l - l29| or in two interacting di- 
atomic molecules [30|, l3l| . Logic gates were driven by 
pulses designed by optimal control or genetic algorithms 
or by using Stimulated Raman Adiabatic Passage (STI- 
RAP) techniques. From the viewpoint of scalability, it 
is challenging to increase the number of encoded qubits 
on a single molecule because the number of cigenstates 
individually addressable cannot grow exponentially. A 
more promising strategy would be to use a network of 
polar molecules holding each a limited number of qubits 
[32 |. Polar molecules are ideal systems for such strat- 
egy for they can interact through the strong anisotropic 
long range dipole-dipole interaction enabling couplings 
between qubits to create entanglement [HI, [33|. Fur- 
thermore, experimental progresses obtained in the for- 
mation of stable heteronuclear alkaline molecules at ul- 
tracold temperatures (~ fj,K) open up the possibilies for 
an individual tight confinement and manipulation of the 
molecules by optical means paving the way to the actual 
realization of a polar molecule based quantum register 
[3^ . Different schemes have been proposed to realize uni- 



versal quantum gates and manipulate qubits encoded on 
molecular levels of ultracold polar molecules. The real- 
ization of a conditional phase gate has been proposed by 
using either a switchable dipole-dipole interaction [34l.l35j 
or a sequence of laser pulses [HI . 

In this article, we present a numerical investigation 
of the possibilities of driving quantum algorithms with 
laser pulses on a register of polar molecules trapped 
in an optical lattice and experiencing a static electric 
field. A step toward scalability is addressed by consid- 
ering schemes that involve intramolecular (implemented 
on a single molecule) and intermolecular (when inter- 
acting molecules are used) quantum gates. Encodings 
that use both the vibrational and rotational modes of 
the molecules are considered as well. We first implement 
the addition of or 1 to a binary digit 6,; and a carry 
in Ci to obtain the sum Si and the carry out c+i [371]. 
This adder requires three qubits and never involves su- 
perposed states. In this sense the computation is weakly 
affected by dephasing. By omitting the phase constraint, 
the conditional population inversions can be realized by 
7r-pulses. Implementing these arithmetic operations on 
molecular systems at this quasi-classical level should al- 
ready be appealing for future applications. This logical 
scheme allows for the concatenation of several arithmetic 
operations by saving the sum Sj in the vibration of one 
molecule for a latter reading and most importantly, by 
transfering the carry out Cj+i to a neighboring molecule 
enabling the next cycle of addition to be carried out with- 
out any intermediary readings. A carefull treatment of 
the phases for quantum gate involving superposed states 
is realized in the context of the Deutsch-Josza algorithm 
[38j applied on a one-bit function. Similarly to the in- 
termolecular gates of the adder, the active qubit and the 
ancillary one are encoded in two neighbouring molecules 
(see also [3(| for a two- molecule implementation) . Pulses 
are then designed by optimal control theory (OCT) with 
an additional constraint for the phase [f| 0] . 

This article is organized as follows: first we give a de- 
tailed description of the polar molecule based quantum 
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register [32| we consider as a support for realistic nu- 
merical simulations and the theoretical framework of our 
treatment. Then we present the conditional population 
inversion driven by it pulses and the realization of a three- 
qubit 0- and 1-adder. Finally, we present the simulation 
for phase correct gates with application to the Deutsch- 
Josza algorithm. 



II. MODEL 

A. Quantum register of polar molecules 

Figure [T] shows a schematic description of the polar 
molecule based quantum register. It consists of a string 
of polar ultracold heteronuclear NaCs dimers cxperienc- 

In the ground X 1 Y 1 + elec- 



ing a static electric field [32 



tronic state, NaCs has a permanent dipole moment of 
4.6 Debye at the equilibrium distance of 7.20 a.u. and is 
among the strongest of all alkaline mixtures [39j making 
it a good candidate for quantum computing. Photoasso- 
ciation from an ultracold Na and Cs mixture can be used 
to obtain NaCs dimers at temperatures cold enough to 
be optically trapped [5c| . Many other ultracold heteronu- 
clear dimers have been obtained using various techniques 
[HI . Experimental techniques for trapping and manipu- 
lating ultracold molecules arc rapidly developping [33J. 
We assume that the molecules can be trapped in the 
lowest translational states of successing sites of a three- 
dimensional optical lattice, with only one molecule per 
site and without tunneling from one site to another. Both 
the string of molecules and the electric field are aligned 
along the Z-axis of a laboratory-fixed frame chosen to 
be the quantization axis. The electric field orientates 
the molecules resulting in a mixing of the rotational lev- 
els. The electric field may vary along the Z-axis to make 
the molecules individually addressable by spectroscopic 
means. 

Molecules are initially stored in a high vibrational 
level of the ground potential. Information can 

be stored in these long-lived levels for they have no per- 
manent dipole moment and the molecules can therefore 
be regarded as isolated from one to another even in the 
presence of the electric field. When a logical operation 
needs to be carried out, the required molecules (depend- 
ing on the number of qubits for the computation) are 
brought to the lowest vibrational levels j35j. As op- 
posed to before, these levels have a large dipole moment 
and the strong dipole-dipole interaction coupling adja- 
cent molecules allow logic gate operations implemented 
on several molecules. One should emphasize that this 
procedure differs from logic gate based on switchable 
dipole-dipole interactions proposed by Kuznetsova et al 
|35j because in our case, the variation of the dipole- 
dipole interaction is not part of the logical operation 
itself. For molecules separated by a distance R of 300 
nm, the dipole-dipole interaction is of the order of a 
few tens of kilo-Hertz. The interaction decreases rapidly 
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FIG. 1: Top figure (a): schematic of a polar molecule based 
quantum register [33]. Bottom figure (b): Permanent dipole 
moments /ij, i=l,2 of the molecules. Their orientations are 
caracterized by the angles 9i and fa. The electric field E is 
parallel to the Z-axis connecting the center of mass of the 
molecules. In the simulation, the wave length of the lattice 
laser is A = 600 nm so the intermolecular distance is 300 
nm, the dipole moment is 4.6 Debye and the dipole-dipole 
interaction is about a few tens of kilo-Hertz 



as the distance between molecules increases due to its 
1/R 3 variation. At 500 nm the interaction magnitude 
drops to a few hundreds of Hertz. Desired driven transi- 
tions arc differentiated from unwanted transitions by the 
dipole-dipole interaction. A small ineraction will require 
a longer pulse duration in order to differentiate the tran- 
sitions. The separation between molecules is dictated by 
the caracteristics of the optical lattice the molecules will 
be loaded in. We choose an intermolecular distance of 
300 nm in our simulations. For NaCs, a wavelength of 
600 nm corresponds to a molecular transition above the 
3s + 7s asymptote and below the inner well of the 1 S + 
potential correlating to the 3s + 7p asymptote [4l| in a 
region were overlap integrals with the lowest vibrational 
levels of the ground 1 £ + are extremely small offering a fa- 
vorable wavelength window for an optical trap with sites 
distant by a distance of 300 nm. Similar features can be 
found for others heteronuclear alkaline molecules |42|. 



B. Hamiltonian 

This work is restricted to the study of gates involv- 
ing two or three qubits encoded in no more than two 
molecules. From now on, we will use capital letters to 
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refer to coordinates in the laboratory frame and mi- 
nuscule letters when we refer to inner coordinates of 
the molecules. We note (ex,ey,ez) the basis vector 
of a Cartesian space-fixed frame (or laboratory Carte- 
sian frame). We note eg,Q = 0, ±1 the basis vec- 
tor for a space-fixed spherical coordinate system. We 
note (e x ,e y ,e z ) the basis vector of a Cartesian molecu- 
lar frame and e p ,p = 0,±1 a molecular spherical frame. 
When the spherical indices are ambiguous, we may use 
the index R and r instead of the index for the laboratoy 
and molecular frame respectively. 

We consider molecules in the ground X 1 Y> + electronic 
potential only. For NaCs, we use the potential calcu- 
lated as explained in [39| . The total Hamiltonian for two 
molecules can be written as the sum of a time indepen- 
dent term Hq and a time dependent pertubation term 
W(t): 

Htot(t) = H + W(t) (1) 

W(t) represents the interaction of the molecules with 
a laser pulse polarized along the ez direction: Ex,(i) = 
E L (t).e z 

W(t) = ]T -/ii.E L (i) = ]T -fi t0 .E L (t) cos 9 t (2) 

i i 

where /ij is the z% = component of the dipole mo- 
ment of molecule i in the ground electronic state. Ho, 
which represents the energy of the molecules coupled by 
the dipole-dipole interaction, is made of several contri- 
butions: 

2 

H = ^(HL + H^ + H^) + V dd (3) 

i=l 

Hg X , H| ra , and Hg are the external, internal, and Stark 
Hamiltonians associated to the molecule i = 1,2 respec- 
tively. V^d is the dipole dipole interaction. 

The external Hamiltonian H l ex describes the motion of 
the molecule i in the trapping optical potential. It can be 
approximated locally by a three dimensional isotropical 
harmonic oscillator with frequency lul and a depth Vq 

M 

HL = ||+V opt (i?;) (4) 

where M is the total mass, ujl and Vq depend on the 
dynamic polarizabilitics of the molecule and on the fre- 
quency and intensity of the trapping laser [H, |43| . Typ- 
ical values for ul are kilohertz. The depth must be of 
the order of tens of kilohertz in order for the molecules 
to be tightly trapped and for the tunneling from one site 
to another to be negligible. 

The inner Hamiltonian H' in describes an isolated di- 
atomic molecule in the Born Oppenheimer approxima- 
tion. It contains the vibrational Hamiltonian H^ ib in the 



ground X 1 Y, + electronic state and the rotational Hamil- 
tonian H!J, ot : 

H in = H vi6 + H rot ( 5 ) 

with 

H ti6 = |^+v(r i ) (6) 

where m is the reduced mass of the molecule vibrating in 
the interatomic potential v(ri), r,; being the interatomic 
coordinate. 

For 1 £ molecules, the rotational Hamiltonian for a 
rigid rotator is given by: 

Kot = BIN 2 (7) 

where N is the total angular momentum and B l v is the 
rotational constant in vibrational state v. 

Ti l s refers to the Stark Hamiltonian of a molecule in a 
static electric field E^(i?j). In our case, ~E l s (Ri) is aligned 
along the Z axis, and varies only along this coordinate. 
The Stark Hamiltonian reduces to: 

& s = -E i (Z i )ii io caB(0 i ). (8) 

The last contribution comes from the long range 
dipole-dipole interaction term V^R.) of two polar 
molecules separated by a distance R = R2 — Ri = 
Rez = Ren in the laboratory frame. By assuming that 
the molecules rotate in a plan, the projection of the rota- 
tional quantum number on the internuclear axis is con- 
served and the dipole-dipole interaction term becomes 

V dd (R, 0102) = - ^ cos 0! cos 9 2 (9) 

27re R 

where the angles 9\ and 62 caracterize the orientation 
of the dipole moments versus the intermolecular axis e#. 

C. Product and coupled basis sets 

For a system of two interacting molecules, we define a 
product basis set \ip) constructed by tensorial product of 
the individual basis of the individual molecules. 

\ip) = \n,v,N,m N ) 1 (g>\n,v,N,m N ) 2 (10) 
= |ni,n 2 ,i>i, N 1 ,m Nl ,V2,N 2 ,m N2 ) (11) 

where the individual basis set |n, v, N, mjv)i is built from 
a tensor product of the cigenbasis of H^, H* ib , and H l rot . 

We make the assumption that ~Vdd and H5 only couple 
rotational levels. However, both terms can create cou- 
plings between states of the optical trap. For a molecule 
in the ground state of the trap, the translational ground 
state wave-function of the potential well is a Gaussian 
of width a = ^JTl/Mlol. For a typical experiment, a is 
of the order of tens of micro-meters. No coupling will 
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occur if H$ is constant over the range of variation of 
a, the coupling integral (nj|Hg|mj) involving the eigen- 
functions of the harmonic trap being zero in this case. 
For the couplings due to the dipole-dipole interaction, a 
detailed analysis can be found in ref |35[. Coupling be- 
tween translational states for typical polar molecules is 
small and results in an energy shift no larger than 1% 
of the unperturbed energy spacing. For typical optical 
trap, where the spacing between translational levels is 
about tens of kilo-Hertz, it gives a negligeable shift of a 
few tens of Herz. the n label can then be omitted in the 
basis function notation. 

Also, we only consider motion in a plane 0, = with a 
field aligned along ez so only = states are involved. 
Finally, the relevant contracted notation \v\, N\, v 2 , N 2 ) 
will label the states of the total product basis. It can be 
further simplified to just |iVi, -ZV2) when different vibra- 
tional states are not considered. 

The states of the coupled basis set that diagonalizes 
the full time independent Hamiltonian Ho are noted 
\<p) = \vi,Ni,V2,N2) (or \(p) = \Ni,N 2 ) when the vi- 
brational manifold is well known) to indicate the product 
basis state Ni, v 2 , A^.to which they are adiabatically 
connected. 
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FIG. 2: Diagram of the lowest eigenvalues of the total time- 
independent Hamiltonian for three different cases: without 
any electric field (a), with electric fields E 1 > E 2 without 
and with dipole-dipole interaction ((b) and (c) respectively). 
The eigenvectors | JVi , JV2) correlate adiabatically to the vec- 
tor I JVi , JV2) with v\ = «a = 0. Only mjv = are considered, 
the 771 at anf m labels are omitted. Transitions in the first 
or second molecule are shown by blue or red arrows respec- 
tively. The subscripts labelling the frequencies indicate the 
rotational state of the neighboring molecule. 



D. Dynamics 

The time dependent Schrodinger equation will be 
solved in the coupled basis set in the interaction rep- 
resentation [ilj]. The evolution of the populations b n (t) 
are given by the set of coupled equations : 

in Jt K(t) = 2 elUJnktw nk(t)h(t) (12) 

k 

where we have introduced the Bohr pulsation u n k = 
(E n —Ej.)/h and the matrix elements W n k (t) of the time- 
dependent interaction in the eigenbasis of Hq . 

W nk (t) = {(p n \W\(p k ) (13) 

The coupled equations are solved using the fourth-order 
Runge-Kutta method (4o| without any rotating wave ap- 
proximation. 



III. INTERMOLECULAR CONTROLLED-NOT 
GATE 

We first examine the possibility of steering elementary 
gates by tt pulses in the microwave and infrared domain 
by focusing on conditional population inversions. The 
control of the phase will be considered in the final section. 

The first example concerns the two-qubit controlled- 
NOT (CNOT) gate which flips the second (target) qubit 
if the first one (control qubit) is equal to 1. 



( a |00)+^|01) + 7|10) + <5|ll»^^ (w) 
(a|00} + j8|01)+7|ll} + <5|10)). 

The gate is driven by a 7r-pulse defined by 

E„(t) = E°sin 2 ^y-j cos(u>i f t) (15) 

with amplitude E® = where r p is the duration of 

the pulse and fii / is the dipole moment for the transition. 
The scheme is similar to some early proposals for condi- 
tional quantum dynamics [!(| [47j . It uses the fact that a 
transition frequency in the second molecule depends on 
the state of the first molecule. 

The qubits are encoded in two neighboring molecules. 
Both molecules are in the ground vibrational states v\ = 
v 2 = so that the states are denoted only by the ro- 
tational states \Ni,Nj). The logical states are encoded 
into the first four states of the coupled basis set which 
correlates with the Ni = 0,1 and N 2 — 0, 1 states (see 
Fig. [2]). In the product basis the CNOT transforma- 
tion should consist in flipping the rotational state of the 
second molecule if and only if the first one is in state 
Ni = 1. In a basis set diagonalizing the Stark Hamil- 
tonian without dipole-dipole interaction, the frequencies 
of the rotational transitions in the second molecule when 
the first one is in the state TVi = (cj ) or Ni = 1 (uji) 
are the same (see Fig[5]). The dipole-dipole interaction 
leads to two very close but different frequencies. This 
allows us to excite only the transition 1 10} —> |11) at a 
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FIG. 3: Evolution of the populations during the CNOT gate 
driven by a 7r-pulse. When the control qubit is (initial states 
1 00) or 1 01) ) , the final state of the qubits remains unchanged 
when the pulse is applied although the population may vary 
during the process (top left and right panels respectively). On 
the contrary, the second qubit is flipped when the initial state 
is 1 10) or 1 1 1 ) leading to final states |11) or |10)(bottom left 
and right panels respectively). 



frequency uj\ (full red arrow in Fig. [5]) and not the tran- 
sition 1 00) — > 1 01) at a frequency £l>o (dashed red arrow). 
When the system is in a computational basis state, a 
7r-pulse can selectively drive the desired population in- 
version. To get a perfect selectivity and avoid the off 
resonance transition the pulse duration must satisfy the 
condition t p > 10/|cl>i — ui \ [48j . 

FigEl shows the evolution of the populations during 
the CNOT gate driven by a 7r-pulse starting from the 
four different states of the computational basis. The 
states corresponding to the logical states 1 00) and |01) 
temporarily leave the computational space but correctly 
return to it with an arbirary phase (the phase constraint 
will be discussed after). From now on, simulations are 
ran with the following fixed parameters: R\2 = 300 
nm, E 1 = 2kV.cm _1 , E 2 =1.5 kV.cm -1 . The dynam- 
ical basis set for the simulation contains the vibrational 
state v = and rotational states up to N — 4 for each 
molecule. The difference of frequencies which must be re- 
solved is 4.002xl0 -6 cm -1 and the pulse duration must 
be t p = 13.2 ^s. The intensity of the 7r-pulse is E® = 7.5 
V.m- 1 . . 



IV. AND 1 ADDER 

The adder of and 1 (0-ADD and 1-ADD) H3 is a 
simplification of the more general full adder of two binary 
digits di and bi and a carry in c,; [49j . Similarly to the 
classical case, the addition of two numbers in binary rep- 
resentation, a = (aia2".<Zi...a n ) and b — (bib2-.-bi...b n ) 
is performed digit by digit starting from the least sig- 
nificant one. For the i th rank, one needs to evaluate a 
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FIG. 4: Logical circuits for the adder (left panel) and 1 
adder (right panel). 



sum Sj and a carry out Cj+i. A cycle of the quantum 
full adder requires four qubits [H, [5(|. However, in the 
special case of a or 1 adders, only three qubits oper- 
ations arc necessary. The C; and bi digits are encoded 
in the first two qubits and the third is 0. One has ini- 
tially \Q1Q2Q3) = \ci,bi,0). The information that a, is 
or 1 repectively is contained in the applied pulse which 
is different for the 0-ADD or 1-ADD gate. The results 
Si and Cj+i are stored in the second and third qubits 
respectively. One has finally \Q1Q2Q3) = \c-i, $i, Cj+i). 
The 0-ADD or 1-ADD drives the unitary transformation 
\ci,bi,0) —t \ci,Si,Ci+i) for a, is or 1 repectively The 
truth tables are given in Table Q] (note that only the first 
four inputs with Q3 = are effectively used, the last four 
ones insure the reversibility of the operations). 

As schematized in FigJU the carry out of the 0-ADD 
gate is Ci. AND. bi and is computed by the three-qubit 
TOFFOLI gate: \x)\y)\z) -> \x)\y)\z ® x.y) where ® is 
the sum modulo 2. The TOFFOLI gate corresponds to 
the AND gate when the third qubit is zero. It is the 
controlled-controlled NOT which flips the third qubit if 
the first two qubits are in state 1. The sum Sj is given by 
the first CNOT gate (noted CNOT1): \x)\y) ->• \x)\y®x) 
The 1-ADD is the continuation of the 0-ADD. For the 1- 
ADD, we can sec from Table Q] that Si is obtained from 
the Si of the 0-ADD by an additional NOT gate which 
simply flips the state (|0) and |1)) of the qubit. It can 
be seen also that Cj+i of the 1-ADD only needs to do the 
sum modulo 2 of the Si and the carry out Cj+i of the 0- 
ADD which are encoded in the second and third qubits. 
As a result this step is a new CNOT (noted CNOT2) 
with control qubit Q2 and target qubit Q3. 

The most promising implementation of the 0-ADD and 
1-ADD consists in encoding the carry Ci and the num- 
ber bi in the rotational and the vibrational structure 
of the first molecule whereas the carry out c,-_|_i is en- 
coded in the rotation of the second molecule. So one 
has \Q1Q2Qz) = \N\,Vy,N 2 )' The sum s, for the i th 
step of the addition will be stored in the vibration of the 
first molecule in replacement of bi and the carry out en- 
coded in the second molecule becomes the carry in for 
the (i + l) th step allowing the addition to be further con- 
tinued. 

When the logical states are mapped on a coupled basis 
set and not on a product basis set, one has to steer the 
transformation in the total computational basis set (a 
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TABLE I: 0-adder and 1-adder truth tables. The information 
cti = and a; = 1 respectively is contained in the pulse driving 
the transformation. Only the first four inputs with Qz = 
are effectively used. Cj+i = 1 © Ci+i. 



three-qubit space here). For a onc-qubit (or two-qubit) 
transformation one has to consider all the transforma- 
tions involving the active qubit(s) in the 3-qubit space. 
For example, the extended operator of the NOT gate on 
Q 2 is 



Uq 2 = E Ql <g> U Q2 <g> E Qs 



(16) 



where EQ k is the unity matrix in the one-qubit space k. 
This leads to four population inversions in place of a sin- 
gle one in a product basis set. Secondly, one should point 
out that whatever the encoding is, unwanted transitions 
can always occur. The length of the pulses are then cho- 
sen to satisfy two opposite conditions. The pulse must be 
short enough to simultaneously drive a maximum number 
of active transitions with a single carrier frequency and it 
must be long enough to avoid these unwanted transitions 
whose frequencies are very close. 



A. Intermolecular TOFFOLI gate 

Conditional dynamics are more demanding for the 
TOFFOLI gate than for the CNOT gate presented in 
the previous section. Although the computional ba- 
sis \Q1Q2Q3) generated by the eight combinations of 
Ni = 0, 1, V\ — 0, 1, N 2 = 0, 1 (with v 2 always equal to 
0) would provide an intuitive mapping between the quan- 
tum numbers and the logical states, it is rather uncon- 
veniant in practice because of the very similar rotational 
constants for the two vibrational levels. With this assig- 
ment, the rotational state of the second molecule must be 
flipped if and only if the first molecule is in a given rota- 
tional and vibrational state. Then, except for N 2 holding 
the carry out, we adopt an assignment with no correspon- 
dance between the quantum numbers and the associated 
qubit state. The eight states and the logical mapping 
are shown in FigJSJ In the present case, the frequency 
of the 1 110) O I 111) TOFFOLI transition (solid arrow 
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FIG. 5: Encoding of the three qubits for the ADDER in the 
states of the coupled basis set \Q1Q2Q3) o \vi, Ni,V2, N2) . 
The active transitions for the TOFFOLI and CNOT gates 
(the latter is called CNOT1 in the text) are shown in full 
arrows and the corresponding unwanted transitions are shown 
by dashed arrows. The zero of energy is chosen at the ground 
state without Stark field and dipole-dipole coupling. 



in FiglS]) is fairly close to the unwanted |000) <-> |001) 
transition. The difference in frequency which must be 
resolved is then 4.0 x 10 -7 cm -1 . 

The population of each computational basis state dur- 
ing the TOFFOLI gate is plotted in the top panel (a) of 
Fig. [6] Each population is weighted by a factor different 
from 1 for clarity (it is not a superposed state). The du- 
ration of the 7r pulse is t p = 253 /is. The carrier frequency 
is in resonance with the active transition 1 1 10) -H- 1 1 1 1 ) , 
i.e. |0, 1, 0, 0) ^ |6, 1, 0, 1) with u= 0.16743 cm" 1 . The 
amplitude of the pulse is Eq = 4.1506 x 10~ 6 kV/cm. 
For all the simulations relative to the adder algorithm, 
the dynamical basis set contains 64 states including the 
i>2 = 2 manifold (v\ = V2 = with J\ and J 2 = 0, 1,2; 
vt =v 2 =2 with Ji = 0, 1, 2, 3, 4). 

The population inversion is very good with a fidelity 
of 0.9999. Note that the tt pulse is not optimized for a 
superposition but is sufficient for arithmetic operations 
for which the system is always in a computational basis 
state. As illustrated in the next section, optimization for 
any superposition can be obtained with optimal control 
theory. 



B. Intramolecular CNOT1 gate 

This gate is intramolecular since both the control Q\ 
and target Q 2 qubits are encoded in the first molecule 
(the rotation and the vibration respectively). The ex- 
tended CNOT1 gate involves the following transitions: 
IIOQ3) <H> IHQ3) for any value of the third qubit Q3. 



7 



•a 



0.2 


noi> 


1100) 


0.15 


"1001) 


HOOO) 


0.1 


1111) 

HI 10) 








1110) 

1111) 


0.05 
°< 


1010) ^iy^>^ 




) 


50 


100 150 


200 




0.2 


1101) 
1100) 








1101) 
1100) 




1001) 










0.15 


HOOO) 










0.1 


1111) 








1110) 


HI 10) 

noil) 








1111) 


0.05 
°( 


ioio) i^SsSrv 




) 


0.05 


0.1 


0.15 





Time (jis) 

FIG. 6: Population evolution for each computational basis 
state during the gates of the ADDER. Upper panel: TOF- 
FOLI gate. Lower panel: CNOT1 gate. Each population is 
weighted by a factor different from 1 for clarity (it is not a 
superposed state) 

These two transitions are shown in full arrow in Fig. [SJ 

wq 3=0 = |100) = |2300) <-> 1 110) = |6l66) (17) 
wq 3=1 = [101) = |2301) O |111) = [0101) (18) 

They correspond to vibrational transitions in the first 
molecule and thus belong to the infrared domain. The 
two transitions being very close (Aw = 8.3 x 10" 7 cm -1 ), 
they can be induced by a single 7r-pulse. Two unwanted 
frequencies are shown by dashed arrows in Fig.©. They 
involve vibrational states of the second molecule. Fortu- 
nately, these transitions differ from the active transitions 
by about 6 x 10" 3 cm -1 . A shorter pulse than for the 
TOFFOLI gate can be used. 

Botton panel (b) of Fig|J5]shows the populations of each 
computational basis state during the CNOT1 gate. The 
pulse duration is t p = 0.194 /is. The carrier frequency 
is fixed by the 1 100) O [ 1 10) transition and is equal to 
wq 3=0 = 196.43 cm" 1 . The amplitude is E = 93.793 
kV/cm. 



C. Intermolecular CNOT2 gate 

This logical operation is IIQ2O) -H> IIQ2I) for any value 
of the second qubit Q^- We then have two active transi- 
tions 

uq 2=0 = 1 100) ee |2300) <-> 1 101) ee |2301) (19) 

UJQ 2 = 1 = 1 110) EE |0106) <-> I 111) EE |6l01) (20) 

Top panel (a) of Figj7] shows the evolution of popu- 
lations corresponding to the intermolecular CNOT gate 
(named CNOT2). A superposition of two microwave ir- 
pulses are used. As for the TOFFOLI gate, we need to 



resolve a frequency difference of 4 x 10 7 cm 1 The pa- 
rameters for the pulses are: t p = 242/is luq 2= o = 0.16734 
cm" 1 , E = 4.1509 x 10" 6 kV/cm and uj q ' 2=1 = 0.16744 
cm" 1 E = 4.1508 x 10" 6 kV/cm. 

D. Intramolecular NOT gate 

The frequencies of the four active transitions of the 
extended NOT gate (\Q1OQ3) <-> \Q1IQ3) for any values 
of Qi and Qz) are 

w Qi= o,q 3 =o = |000) ee I2266) <-> |000) EE |6666) (21) 

wg 1= i,Q 3 =o = |100) = |2306) O |110) ee [0100) (22) 

w 0i =o,q 3 =i = |001) ee |2201) o |011) ee [0001) (23) 

wqi=i,<Js=i = |101) = |2301) |111) = 10101) (24) 

Bottom panel (b) of FigJT] shows the population evolu- 
tions for the final NOT gate. 

The frequencies wq 1= o,q 3 =o and uiq 1= i,q 3= o on one 
side and wq 1= o,q 3 =i and Wq i =i,q 3 =i on the other side 
are sufficiently close to be driven by a single 7r-pulsc 
(Aw = 4 x lb -6 cm" 1 and Aw = 8 x 10" 7 cm" 1 re- 
spectively). The two pulses have a duration adapted 
to avoid unwanted transitions (Aw = 6 x 10" 3 cm" 1 ), 
T p = 0.726/is. The frequencies are wq 1= o,q 3 =o = 196.28 
cm" 1 and wq 1= i.q 3=0 = 196.43 cm" 1 . The amplitudes 
are E = 25.0115 kv/cm and E = 14.4558 kv/cm. 

E. Intramolecular initialization 

Finally, one has to discuss the initialization of an ad- 
dition cycle. The digit bi must be encoded in the first 
molecule for any value of the carry in Cj which is un- 
known. With the assignment used here (see Fig[5]), the 
digit is encoded in the vibration of the first molecule. If 
bi = 1, the molecule must be in the manifold V\ = 0. This 
means that the system is directly ready since the algo- 
rithm generates the carry out in this ground vibrational 
state. If bi = 0, the molecule must be in the manifold 
Vi = 2. Two transitions must then be driven: 

[010) EE |6666) -> 1 000) EE |2206) (25) 
[110) ee |6l66) -> 1 100) ee |2300) (26) 

They correspond to two transitions with frequencies 
wq 1= o,q 3 =o an d wq 1= i j q 3= o of the NOT gate realized 
above and the NOT pulse can then be used to initialize 
the addition cycle. 

F. Concatenation 

Fig. [5] illustrates the full 1-ADDER operation for the 
example b- t = 1 and a = 1. The initial logical state is 
1 110). The network TOFFOLI (intermolecular), CNOTl 
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FIG. 7: Population evolution of each computational basis 
state during the supplementary gates of the 1 ADDER. Upper 
panel: CNOT2 gate. Lower panel: NOT gate. Each popula- 
tion is weighted by a factor different from 1 for clarity (it is 
not a superposed state) 
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FIG. 8: Population evolution of the computational basis state 
corresponding to the logical state 1 110) (pi = 1 and d — 1 dur- 
ing the 1-ADDER gate (TOFFOLI-CNOT1-CNOT2-NOT 
gates) to give the final logical state | 111) corresponding to 
Si — 1 and Ci+i = 1. 



(intramolecular), CNOT2 (intermolecular) and NOT (in- 
tramolecular) gates operating in the microwave or infra 
red domain drives the system towards the final | 111) log- 
ical state corresponding to s, = 1 and Cj+i = 1. 



G. Discussion 

Any perturbation to the energy level structure com- 
ing from neighboring molecules, hyperfine interactions, 
coupling with the translational states in the lattice will 
drastically affect the high fidelity logical operations the 
pulses were designed to drive. 

We chose the electric field to be very different from one 



molecule to the other (2kV/cm and 1.5kV/cm in our sim- 
ulations). The reason for this choice was motivated by 
the fact that we have intramolecular gates implemented 
on one single molecule and well differentiated molecules 
allows for fast intramolecular gates. Having such a large 
variation of the electric field may generate couplings be- 
tween translational states of the trapping potential be- 
cause the requirement of a constant field over the range 
of variation of a might be experimentally challenging to 
achieve. A much smoother gradient might prove to be 
more suitable to this respect. We tried simulations with 
different gradients and it was always possible to drive 
algorithm although intramolecular gate times had to be 
longer when the gradient was small (E 1 = 2 kV/cm and 
E 2 = 1.9kV/cm). On the contrary, the time scale for 
intermolecular gate remains the same being mainly fixed 
by the dipole-dipole interaction. 

We assumed that the electric field was orientated along 
the intermolecular axis. Note that results will remain 
valid for other field orientations. Only the expression of 
the interaction will differ from what we presented in this 
model, but the order of magnitude of the interaction will 
be about the same. 

The effect of neighboring molecules interacting on the 
two-molecule systems has been evaluated. We have 
analyzed the energy shifts on the levels of the active 
molecules j and j + 1 by adding two molecules (at po- 
sition j — 1 and j + 2) and using the same gradient for 
the Stark field as for our simulations. The basis set is 
composed of manifolds v = with N = 0, 1, 2 and v = 2 
with N = 0, 1, 2, 3, 4 for each molecule. The shifts indro- 
duced are of the order of a few kilo-Hertz, the effect is 
only a shift in frequency but the general structure of the 
schemes is not affected. Inter and intramolecular gates 
could still be driven with similar pulses in a larger net- 
work. 

The hyperfine structure of rovibrational molecules 
can complicate the manipulation of molecules with mi- 
crowave pulses because unwanted transitions can occur 
[SlL [52j . In the present model, the hyperfine structure 
was not taken into account. Within the rotational level 
N = manifold, the most important term is the scalar 
hyperfine coupling proportional to ii.i2- where ii and 12 
are the nuclear spins of the atoms of the molecule. The 
splitting between hyperfine levels is typically of the or- 
der of a kilo-Hertz. For N different than 0, the most im- 
portant contribution to the hyperfine Hamiltonian comes 
from the nuclear quadrupolar interaction which gives a 
splitting of the order of a few hundreds of kilo-Herts. In 
our case, = 3/2, and ic s = 7/2. iV=0 splits into 
32 states, and N—l splits into 96 levels. For a static 
electric field of 2kV/cm, more that 4 rotational levels 
per vibrational manifold are needed to described accu- 
rately the energy structure. Taking into account all the 
hyperfine levels would complicate the theoretical simula- 
tions since we had to manipulate the tensor product of 
the individual basis. However, the hyperfine structure 
won't change the conclusion of this work. Simulations 
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of the logical operations showed two different caracter- 
istic time scales for the logic gates. Fast gates can be 
driven with pulses no longer than a few tens of ns. These 
pulses are broad enough that all hyperfine transitions 
within one rotational level will be excited equally, and 
the unresolved hyperfine structure can be neglected. On 
the other hand, slower gates such as the TOFFOLI of 
the adder algorithm are driven by much narrower pulses, 
and the splitting between hyperfine levels, larger than 
the dipole-dipole interaction splitting, will be well re- 
solved and only one hyperfine level will be active. We 
note as well that as the static electric field increases ori- 
entating the molecules along the electric field, mjy be- 
comes a fairly good quantum number, and the matrix 
of the transition dipole moments between hyperfine lev- 
els of different rotational manifold show that for electric 
field typically stronger than lkV/cm, only a few transi- 
tions are allowed reducing considerably the possibilities 
for unwanted transitions. 

Other terms that we have not considered are the second 
order Stark effects resulting from the static electric and 
the optical trapping potential. The static second order 
Stark effect is of the order of a few kilo-Hertz for the 
diagonal terms. The effect due to the laser field can be 
more important due to the wavelength dependence of the 
molecular polarizabilitics. For a laser intensity typically 
of the order of lkW/cm 2 , and a molecular polarizabilty of 
a few thousands of atomic units, the second order Stark 
effect is of the order of one hundred of kilo- Hertz for the 
diagonal terms. It is not strong enough to affect our 
results. 

The numerical simulations suggest that ultra cold 
trapped polar molecules are promising for the concate- 
nation of several gates with a high fidelity. First, we can 



compare the efficiency of a couple of diatomic systems 
with a tetraatomic molecules to implement arithmetic 
operations. It was not easy to encode the four qubits 
of the full adder in the two interacting dimers. Some 
hundreds of microseconds to add 1 is finally very longer 
than the timescale of the full addition in a polyatomic 
entity (some tens of ps). However the main point is 
the possibility to avoid intermediary reading out and 
re-encoding of the carry out. Here the system is directly 
ready for the next cycle involving the next molecule. 
This is a step towards scalability but at the price of long 
pulse duration. The other advantage is the simplicity 
and robustness of the scheme based on 7r pulses. 



V. PHASE CORRECT GATES 

The relative phase of the gate transitions optimized 
among the states of the computational basis set can re- 
duce the accuracy of the gate when it is applied to an 
arbitrary superposed state, ir pulses are not sufficient 
to give phase-correct quantum gates. Then correct gate 
pulses can be determined by the optimal control theory 
generalized to the multi-target case with a phase con- 
straint fiHIHi. 



A. Multi-target optimal control 

The optimal field maximizes the objective functional 
J with the constraint that the Schrodinger equation is 
satisfied at any time [EiLHEl- The functional reads: 



J 



J = £|<^(*/M 



/)| 2 



25R 



W(t)\dt + -H|Vr(i)> 



"t f 

v I E 2 (t)dt 
'o 



(27) 



where a is a positive penalty factor chosen to weight 
the importance of the laser fluence. For a N-qubit gate, 
Z = 2 N + 1 where 2 N is the number of input-output 
transitions in the gate transformations and the supple- 
mentary equation is the phase constraint. The ip™{t) are 
the wave packets which are propagated forwards in time 
with the initial conditions ijj^(t = 0) = (j>f, n = 1,...Z. 
The Lagrange multipliers are propagated back- 

wards in time with the final conditions ip'f(t = tf) = (ftj, 
n = 1,...Z. The supplementary transfer which imposes 
the phase correction is a sum over all the transitions of 



the gate [I? 



fc=i 



5> 



(28) 



and the single phase </> can take any value between and 
2ir. The universal gate field is a sum of Z contributions 



Ej(t) = -(s(t)/ha)S 



£<V>?(*)MC(*)> 



n=l 



(29) 



where j denotes the polarization direction of the electric 
field. The fidelity well adapted to take into account the 
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FIG. 9: Logical circuit for the Deutsch-Jozsa algorithm 



phases is given by: 
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(30) 



B. Deutsch-Josza algorithm 

The Deutsch-Josza algorithm for a one-qubit function 
will serve us to illustrate the realization of phase correct 
gates 

The Deutsch-Josza's problem illustrates the speedup 
of quantum computing by taking advantages of super- 
posed states. It has been the subject of several theoretical 
studies [l3, [H, E3, S. 24 L 3(| , and has been experimen- 
tally implemented (^561- 58^ The principles can be sum- 
marized as follows: we suppose that a function applied 
on one qubit |a;) can either change its value (balanced 
function) or leave it unchanged (constant function). The 
problem is to determine whether the function is balanced 
or constant by a single call to the function and one mea- 
sure. The logical circuit for the Deutsch-Jozsa algorithm 
is sketched Fig. [5] The algorithm requires an auxiliary 
qubit \y). After a NOT gate on \y), the two qubits are 
put in a superposed state by one Hadamard gate on each 
qubit, The latter superposes the qubit states according 
to 



|0) 

ID 



^ (VV2)(|0) 

HAD 



|1» 

(1/V2)(|0) - |1>) 



(31) 
(32) 



The call function is then implemented by applying the 
transformation [//: 



U f :\x)\y)->\x)\y(Bf(x)) 



(33) 



A final Hadamard gate is carried out on |a;). When the 
initial state is 1 00) , the final state of \x) is ±|0) for a 
constant function, and ±|1) for a balanced function. The 
nature of the function is therefore determined by only one 
query. 

We simulate the example with the balanced function 
/(0) = 1 and /(l) = 0. In this case U f is the CNOT gate 
with control qubit |a;). Qubits \x) and \y) are encoded 
in the first two rotational states of two neighbouring 
molecules in the vibrational state v = 0: \x)\y) <-> \N1N2} 
with vi = V2 — 0. (see Fig. The basis set for the 
simulation contains the vibrational state v — and ro- 
tational states up to N = 4 for each molecule. Both 
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FIG. 10: Evolution of the populations during the gates involv- 
ing HADAMARD gates in the Deutsch-Josza algorithm. Top 
panel: NOT-HADHAD step (NOT gate on \y) and the two 
HAD gates on |a;) and \y). Bottom panel: HAD gate on \x). 
The initial state is a superposed state with equal weights 0.25 
on the computational basis set states 1 00) , |01), 10), |11). 



Hadamard steps involve intramolecular transitions and 
the frequencies are mainly fixed by the Stark levels so 
that the pulse duration is shorter than for the intermolcc- 
ular CNOT gate which is the bottleneck of the full algo- 
rithm. 



C. NOT-HADHAD step 

Optimal control allows us to optimize a global gate 
driving the resulting transformation of the first NOT gate 
on \y) and the two HAD gates on both qubits \x) and 
\y) in a single shot as suggested in previous works (25l - 
l28j . This provides an interesting speedup. This NOT- 
HADHAD transformation in the coupled basis set is 

|65> -> 2" 1 (|66) - 1 01> + 1 10) - |11)) (34) 

[01 > -> 2- 1 (|66) + |6l) + 1 10) + |11)) (35) 

|I6> -> 2- 1 (|6o) + |6l) - 1 10) - |11)) (36) 

|11 > -> 2- 1 (|00) - |01> - 1 10> + |11)) (37) 

and the additional phase equation is 

2~ 1 (|60) + |01> ) + [10) + |11)) -> |00) (38) 

The NOT-HADHAD gate has been implemented by op- 
timal control without guess field and a pulse duration 
t p = 63 ns. A fidelity of 0.99999 is reached after 400 iter- 
ations. The phases are optimized within 3 x 10 _3 7r. The 
evolution of the populations starting from a superposed 
state with equal weights 0.25 on the computational basis 
set states [00), |01), |10), |11) is shown in top panel of 
Fig 1101 The corresponding optimal field and its Fourier 
transform are given in Fig. (|lll) . 
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FIG. 11: Amplitude of the optimal field for the NOT- 
HADHAD gate (left top panel) and for the HAD gate (right 
top panel). Bottom panel; the corresponding Fourier trans- 
forms. 



D. Phase correct intermolecular CNOT gate 

The supplementary transformation to ensure that the 
phases of the final states of each transition are equal is 
here 

i(|oo> + |oi) + |io> + |ii»-> 

i(|00)e^ 1 + |01)e^ 2 + |10)e^ 3 + |ll)e^ 4 ) (39) 



with 



(40) 



where the phase tp can take any value between and 
2ir. The ir pulse is a very good trial field since the conver- 
gence is fast with a fidelity index of 0.99943 in 130 itera- 
tions. Fig |12l compares the population evolution starting 
from a superposed state driven by the ir pulse (see Figj3]) 
or by the optimal control field. The fidelity for the trans- 
formations of the states |00), 1 01) decreases with the ir 
pulse. Owing to the phase constraint, the fidelity is bet- 
ter with optimal control. The phases obtained for the 
fifth transformation (|39|l are ipi = 0.2027T, tp 2 = 0.1997T, 
(ps = 0.2007T, and tp 4 = 0.2027T. The optimal control has 
found field amplitudes of the same order of magnitude 
as for the 7r-pulse. The main frequency is the carrier fre- 
quency of the 7T pulse and a lot of very small frequencies 
without direct signification from the model. 



E. HADAMARD gate on \x) 

One has to consider the extended transformation in a 
two-qubit space here. The HAD gate (last step) on \x) 



FIG. 12: Evolution of the populations during the phase 
CNOT gate starting from a superposition of states |00 >, 
j 1 >,|10 >, 1 11 >. Top panel: gate driven by the 7r-pulse 
calculated in Fig. The pulse optimized for a pure state 
does not lead to a high fidelity operation when applied to a 
superposition of states; bottom panel: high fidelity gate opti- 
mized by optimal control. We checked that high fidelity gates 
were obtained when the pulse were applied to other initial 
superpositions. 
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FIG. 13: Amplitude of the optimal field for the CNOT gate 
(top panel) and its Fourier transform (bottom panel). 



becomes in the coupled basis set: 



|00) 
|01) 

|io> 

111) 



2- 1 / 2 (|6o) + |6i)) 
2- 1 / 2 (|6i) + |n» 
2" 1 / 2 (|06) - | io> ) 
2- 1 / 2 (|6i)-|ii)) 



(41) 
(42) 
(43) 
(44) 



with the additional equation for the phase constraint. 

(|00) + |01> + 1 10) + |11)) ->■ 2~ 1 / 2 (|06) + [01» (45) 

The evolution of the populations starting from a super- 
posed state with equal weights 0.25 on the computational 
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basis set states |00), |01), |10), |11) is shown in bottom 
panel of Fig [TDJ The corresponding optimal field and 
its Fourier transform are given in Fig fTl] As before, a 
fidelity of 0.99998 is obtained after 400 iterations. The 
phases are optimized within 10 _2 7r. 

Final optimal fields driving the desired operations can 
have a complicated envelop and differ drastically from the 
initial 7r-pulse raising the question of the robutness and 
actual realisation of such pulses. For schemes involving 
only rotational levels, the Fourier transform of the pulses 
show that the spectrum contains frequencies in the mi- 
crowave region only. In the case of the HADAMARD 
gates (see FiglTTj). frequencies are within 0.1 cnr 1 (3 
GHz) and 0.2 cm -1 (6 GHz)). For this range of frequen- 
cies, electronic equipments can generate any arbitrary 
pulses and therefore obtaining such optimized pulses is 
experimentally feasible. Because several rotational levels 
are taken into account in the model (up to N=A in the 
present study) , optimized pulses may contain frequencies 
that will be more challenging to obtain. In the case of 
the optimized CNOT gate (see Fig fT5|) . the pulse con- 
tains frequencies up to 0.75 cm" 1 (20GHz). Filters can 
be used during the optimization procedure to eliminate 
out of range frequencies. If different vibrational levels are 
used. Pulses will be in the infrared region where pulse 
shaping techniques are less developed. 

VI. CONCLUSION 

Manipulating trapped ultracold molecules by laser 
fields offers great potentialities toward scalable quantum 



computing. Logical operations can be performed by 
splitting processes into intramolecular global gates 
and intermolecular gates. The latter are very crucial 
for scalable operations for they enable the transfer of 
informations from one molecule to another through the 
dipole-dipole interaction. Using the states of the coupled 
basis set as logical states complicates the realization of 
one-qubit or two-qubit gates when the number of qubits 
involved in the logical operations increases. This can 
be overcome by using schemes that involve switchable 
interactions. Storing molecules in rovibrational levels 
with vanishing average dipole moment and transfering 
them when needed to levels with strong dipole moments 
by adiabatic passage techniques could be a possibility. 
This operation would not take any longer than a 
few microsecond and would be equivalent in terms of 
duration to an additional intermolecular gate. The 
implementation of schemes that would include switch- 
able interactions would further extend the potentialities 
of polar molecules for scalable quantum information. 
This remains open and can stimulate new researchs in 
molecular design, storage and communications among 
entities. 
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